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, Abstract 

We investigate some exact static cylindrically symmetric solutions for a perfect fluid 
O ' in the metric f{R) theory of gravity. For this purpose, three different families of so- 

lutions are explored. We evaluate energy density, pressure, Ricci scalar and functional 



^ ' form of f{R)- It is interesting to mention here that two new exact solutions are found 

from the last approach, one is in particular form and the other is in the general form. 
The general form gives a complete description of a cylindrical star in f{R) gravity. 
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^ I 

■ 1 Introduction 

O 

' Latest data from CMBR and supernovae surveys show that the energy composition of the 

universe is as follows: 4% ordinary brayonic matter, 20% dark matter and 76% dark energy 
Dark energy seems to be the most favorite way to explain the accelerated expanding 
^ ' universe as indicated by recent observations. However, the exact nature of dark energy is 

^ ■ still unknown. It is assumed to have a large negative pressure. There are some problems 

in astrophysics and cosmology like the issues of dark energy and dark matter where we 
experience severe theoretical difficulties. 

Many proposals are available in literature [5]-[9] to generalize the theory of general 
relativity. The f{R) theory of gravity is the simplest generalization such that the Ricci 
scalar in Einstein-Hilbert action is replaced by an arbitrary function of the Ricci scalar. 
This generalization leads to complicated field equations as well as of higher order. Due 
to higher order derivative dependence, we expect to get more exact solutions than general 
relativity. This theory explains the accelerated expansion of the early universe as well as 
the current accelerated expansion of the universe. 

Recent literature indicates keen interest to investigate exact solutions in f{R) gravity. 
Multamaki and Vilja |10j discussed static spherically symmetric vacuum solutions. The 
same authors jllj extended it for perfect fiuid and found that pressure and density cannot 
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be determined uniquely unless some conditions are imposed. Carames and de Mello |12j 
explored static spherically symmetric vacuum solutions in higher dimensions. Capozziello 
et al. p!3] found static spherically symmetric vacuum solutions via Noether symmetry 
approach. Hollenstein and Lobo [U] analyzed exact static spherically symmetric solutions 
by coupling f{R) gravity with nonlinear electrodynamics. Sharif and Shamir [T2] explored 
static plane symmetric vacuum solutions. The same authors [16] found exact solutions for 
Bianchi types I and V spacetimes in vacuum. Recently, Sharif and Kausar [17] discussed 
static spherically symmetric dust solutions. The same authors [18j also studied non-vacuum 
solutions of Bianchi V/q universe. In a recent paper, Shojai and Shojai [19] found some 
static spherically symmetric solutions for perfect fluid satisfying the physical acceptability 
criteria. 

The study of exact solutions in f{R) gravity is mostly restricted to spherical symmetry 
either vacuum or non-vacuum. One can compare the results with solar system observations 
based on the Schwarzschild solution. The physical acceptability of these solutions can also 
be discussed. Azadi et al. [2^ investigated static cylindrically symmetric vacuum solutions. 
Momeni and Gholizade [21] found static cylindrically symmetric vacuum solutions with 
constant scalar curvature applicable to outside of a string. It would be interesting to extend 
the study of exact solutions for cylindrically symmetric non- vacuum solutions. Delgaty and 
Lake [22] gave the physical acceptability criteria for an exact solution in GR as follows: 

• isotropy of pressure; 

• regularity at the origin; 

• positive definiteness of the energy density and pressure at the origin; 

• vanishing of pressure at some finite radius; 

• the monotonic decrease of the energy density and pressure with increasing radius; 

• subluminal sound speed, i.e., (u^^ = ^ < 1). 

For f{R) gravity, only the first and last two conditions are accounted. 

In this paper, we construct some static cylindrically symmetric interior solutions in f{R) 
gravity. We consider three types of assumptions depending upon model parameters to ex- 
plore these solutions in tabular form and discuss them by using the physical acceptability 
criteria. In particular, we are interested in physically acceptable solutions of a cylindrical 
star. We discuss the dependence of f{R) functions on the scalar curvature for these solu- 
tions. The outline of the paper is as follows: In section 2, we give a brief overview of the 
metric f{R) gravity and formulate the field equations of cylindrically symmetric spacetime. 
Section 3 provides solutions of the field equations by using three types of assumptions for 
particular values of the parameters, given in the form of tables. In the last section, we 
discuss and summarize the results. 
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2 Field Equations in f{R) Gravity for Cylindrically Symmet- 
ric Spacetime 

The modified form of the Einstein-Hilbert action in f{R) gravity is given by [23] 

S = I d^xV^[f{R) + KCm], (1) 

where g is the trace of the metric tensor g^u, f{R) is a generic function of the Ricci scalar 
R, K is the couphng constant in gravitational units and Cm is the standard matter La- 
grangian. Variation of the action with respect to the metric tensor leads to following fourth 
order partial differential equations 

FiR)R^, - ]^f{R)gt,y - V,,VuFiR) + g^,yUF{R) = kT^,, (2) 

where T^jy is the energy-momentum tensor of matter, F{R) = df{R)/dR, □ = V^V^ with 
representing the covariant derivative. Taking trace of the above equation, we obtain 

F{R)R - 2f{R) + 3DF{R) = kT. (3) 

Here R and T are related differentially and not algebraically. This indicates that the field 
equations of f{R) gravity will admit a larger variety of solutions as compared to GR. Notice 
that R = implies T = in GR whereas in f{R) theory, this does not hold. The Ricci 
scalar curvature function f{R) can be expressed in terms of its derivatives as 

^ -'^T + F{R)R + 3aF{R) ^ 

Substituting this value of f{R) in Eq.(l2|), we obtain 

RapF{R) - ^[RF{R) - UF{R)]g^p - V^VpF{R) = k[To,^ - \iP - ^P)9ap]- (5) 

The obtained equation is independent of f{R) and is used to calculate the independent field 
equations of a system. 

The line element of static cylindrically symmetric spacetime is given by [23] 

ds"^ = A{r)dt^ - B{r)dr^ - r^{de^ + a^dz^), (6) 

where A, B are the metric coefficients and a has dimensions of 1/r. The energy-momentum 
tensor for perfect fluid is 

= {p + p)ui^Uu - pg/Mu, (7) 

here = is the four-velocity in co-moving coordinates, p is the density and p is the 
pressure of the fluid. The corresponding Ricci scalar has the form 



AA' X' (A! 2A\ 2A 

A + + + 



r X \ 2 r j r 



(8) 
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where X = AB and prime denotes derivative with respect to the radial coordinate r. Using 
Eq.©, we get the following two independent equations 

X' , X' 2rKX , 
-2rF" + r—F' + 2F—-^-{p + p) = 0, (9) 



F' X'\ f 2A\ 2A 2kX 

The conservation equation, T^.^ = 0, leads to 

A' 2p' 



a"+{---](^'- — ]-z:t-^(p + p) = 0- (10) 



(11) 



A p + p 



Thus we have a system of three non-linear differential equations with five unknown functions, 
namely, F{r), p{r), p{r), A{r) and B{r). 



3 Solutions of the Field Equations 

This section is devoted to discuss the solutions of the field equations. The Ricci scalar 
depends on A{r) and B(r), so it would be much difficult to solve the field equations involving 
f{R). Consequently, analytical solutions of Eqs.(f9|)- pT]) could not be found in closed form. 
Here, we use Tolman's method |25j which constructs a variety of explicit solutions of the 
Einstein field equations. Tolman introduced a relation for the metric components depending 
on radial coordinates to solve the set of equations for perfect fluid. This procedure yields 
pressure and density such that they may or may not be physically acceptable. We explore 
cylindrically symmetric exact solutions in terms of energy density, pressure and functions 
of R. In order to solve the field equations, we use three types of assumptions [TOl [H] for X 
and F given below: 

. X = Xo, F = For-, 

. X = Xor™, F = Fo, 

• X = Xor™, F = For"-, 

where Xq, Fq, n and m are arbitrary constants. For any choice of m and n, a number of 
interesting solutions for A, p, p and R are obtained along with the corresponding function 

3.1 Type I: X = Xq, F = Fqt'' 

Inserting X = Xq, F = Fgr" in Eqs.(l9|) and (jlOp and combining the resulting equations, 
we obtain 

, „ n 2r? — An — 2 , , , 

A" + -A' + -2 A = 0. (12) 

Its general solution is 

A = Aor™, m=^(l-n± ^9 + Un - Tn^) , (13) 
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where is a constant. Using this value in Eqs.(j9|)- ()lip . it follows that 

FoAQmn{n — 1) 



P = Po +Pir 
where po is another constant. Further 

The Ricci scalar takes the form 



m+n—2 



Pi 



Pi 



i?1 



2KXo(n + m-2)' 

FQAQn{n - l)(2n + 3m - 4) 
2KXo(n + m-2) 



+ 3m + 2). 



(14) 

(15) 
(16) 



We can write some exact solutions for particular values of m and n given in Table 1. The 
last two columns have functional form of f{R) and speed of sound, respectively. We have 
already assumed the isotropy of pressure. This type does not provide any solution which 
satisfies the physically acceptable criteria completely as the solutions I, III, V, VI and VII 
indicate squared sound velocity as negative. The solutions II and IV do not have monotonic 
decrease in energy density and pressure with increasing radius. 

Table 1: Type I solutions with X = Xq, F = For"^. 



No. 


m 


n 


A{r) 


p{r) 


p{r) 


R{r) 


f{R) 




I 


-3 


1 


Aor 


Po + ^ 




dAo 
Xnr 




^3 


II 


-2 





Aor' 


Pa 


-PQ 


l-2Ao 
Xn 


any regular 


< 1 


III 


-2 


1 




Pq 


-Pa 





any regular 


-1 


IV 





1 


Aor 


Po 


-Po 


(-iAo 
Xnr 


any regular 


< 1 


V 





2 


Aor' 


Pq 


-Pq 


l-2Ao 
Xn 


any regular 


-1 


VI 


2.5 


1 


Aor 






6Aa 


n Fo(3i?-i/3_5^)i/2 


1 

4 




^0 ^ 2kXo 


Xnr-* 


Jo 6 


VII 


2.5 


2 


Ao 
r 


IhAoFo 


PO ^ iKXoV^ 





any regular 


1 
2 



3.2 Type II: X = Xor™, F = Fq 

Now replacing these values in Eqs.Q and (jlOp . it follows that 

A" 



^ A' 2 + ^ 

2r r'^ 



which leads to the following solution 



A = Aor"^, n=j(2 + m± y/ m'^ + 20m + 36 



Using this solution, we can write pressure, density and Ricci scalar as 

n-m-2 FoAomn 



p 
p 

R 



Po+Pir ' '■- Pi 

n—m—2 



Po + Pir 

n—m—2 



Pi 



Rir 



2nXo{n — m — 2) ' 
FoAom{3n - 2m - 4) 
2KXo{n — m — 2) ' 

Ri = -r^(n^ + 3n — 2m — nm + 2). 
5 



(17) 

(18) 

(19) 
(20) 
(21) 



For the same choices of m and n taken for type I solutions, we construct a Table 2 for 
this type of solutions given below. Following the physical acceptability criteria, we see that 
solutions II and V give negative squared sound velocity, while it becomes 1 for solution IV. 
Solutions I and III do not show decrement in energy density and pressure by increasing the 
value of radius. Only solution VI satisfies complete criteria but the Ricci scalar turns out to 
be radial dependent. Hence, we could not obtain any acceptable solution for the constant 
curvature condition. 



Table 2: Type II solutions with X = Xor"^, F = Fq 



No. 


m 


n 


Mr) 


P{r) 


p(r) 


i?(r) 


f{R) 


y 2 _ dp 

'^s ~ dp 


I 


-2 





Ao 


Po 


-Po 


Ma 

Xn 




< 1 


II 


-2 


1 


Aor 


Po 


1 Aol-hr 
^ i^Xn 


~P0 


kXc 


V2Aar 
X" , 


fo + 


1 

3 


III 


-3 


1 


Aor 


Po - 




-pa - 




15Aor^ 




< 1 


4K.Xn 


4K.Xn 


Xn 




IV 





1 


Aor 


Po 


-Po 


Ma 

Xnr 


any regular 


1 


V 





2 


Aor^ 


Po 


-po 


l-2Ao 


any regular 


-1 


VI 


2.5 


2 


Aor' 


Po ^ 


AoFo 


-Po - 




2Ao 




< 1 


'kKXar-'l'^ 




Xar-'/^ 


J" 20 



3.3 Type III: X = Xor", F = For"" 

This type of assumption deals with radial dependent expressions both for X and F, i.e., 

X = Xor"" F = For"" 
Replacing these values in Eqs.(l9|) and p^ . it follows that 



, „ 2n — m , 2'n? — An — m — 2 — nm , 
A" - — A' 2 ^ = 0' 

which has the following general solution 

A = Air'i + A2r'2, 

where Ai, A2 are constants and /i, I2 are given as 
1 



(22) 



(23) 



li = — — 2n + \/32 + 16m + m? + 64n + 12mn — 2872^ 
^2 = - (^rn — 2n — \/ 32 + 16m + + 64n + I2mn — 28n2 



Inserting the above values in Eqs.(j9])- (fTT|l . the corresponding p, p and R are obtained as 
with 



FoAi/i(2n^ - 2n - 2m - mn) Fo^2^2(2n2 - 2n - 2m - mn) 

Pi = . ^ n . > P-2 = 



AkXq{Ii + n - m - 2) 



4:kXo{12 + n - m - 2) 



Also, 

P=-P0 + pir'i+"-"^-2 ^ ^^^l,+n-m-2^ (25) 

where 

_ FoAi (3/i + 2n - 2m - 4) (271^ - 2n - 2m - mn) 

~ 4:kXo{Ii + n - m - 2) ' 

_ F0A2 (3/2 + 2n - 2m - 4) (2n2 _ 2n - 2m - mn) 
~ 4«;Xo(l2 + w-TO-2) ■ 

The corresponding Ricci scalar is 

R = i2ir'i-"»-2 + i?2r'2-"^-2, (26) 

where 

Ri = 7^(2/? + (6 - m)h + 4m), i?2 = 4^(2^1 + (6 - m)Z2 + 4m). 

We formulate solutions corresponding to Zi and I2 in two different tables. Firstly, we calcu- 
late solutions for Zi given in Table 3. 



Table 3: Type III solutions with X = Xqt'^, F = For'^ corresponding to h. 



No. 


m 


n 


A{r) 


p{r) 


p{r) 


R{r) 


m) 


-s' = Z 


I 


-12 


-1 


Al 

r 


^0 2kXo 


■f^O 2kXo 


-32Air" 
^0 


/o - aFoi?8/9 
where 

" W9 V Xn ' 


1 

8 


II 


-10 


-1 


Ai 


PO 5kXo 


PO 5kXo 


-32/lir" 
Xo 


fo - a6Foi?5/6 
where 


1 

4 


III 


-2 





Ai 


Po 


-Po 


-4Ai 


any 


< 1 


IV 


1 


3 


Ai 


Po 


~Po 


2Ai 
X„r:' 


any regular 


< 1 


V 


4 


1 






6A,F„r-^ 
P° kXo 


Xor^ 


where = 

V -^0 


2 


VI 


6 


2 


Air 5 




15A,F„r-' 
PO kXo 


37Ai 
Xor-> 


fo + abFo^ 
6= (4;;)2/3 


< 1 


VII 


8 


5 






PO kXo 


26Ai 


/o + a6Foi?i/6 
where 

n — 51 

^= (Al)5/6 
V An ' 


2 



Since we have already assumed isotropic pressure, so this condition automatically holds 
for all the solutions. We see that solutions III and IV have subluminal speed but energy 
density and pressure do not decrease with the increase in radius. The solutions I and II show 
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decrease in density and pressure as radius increases but these give negative squared sound 
velocity. For solutions V and VII, the energy density and isotropic pressure decrease when 
the radius increases and > 1. The solution VI fulfills the criteria of physical acceptability 
properly. Hence this is the only acceptable solution obtained by using this assumption. 

Table 4: Type III solutions with X = Xqv"^, F = For'" corresponding to l2- 



No. 


m 


n 


A{r) 


p{r) 


p{r) 


Rir) 


fiR) 




I 


-2 





M 

r 


Po + ^^ 


^2 . 


-7A2 

Xnr 


fo - j^R 


< 1 


II 


4 


1 






PO ^ 8kXo 


I4A2 
Xor'' 


where 

n - 91 

" — 226/9x7«/9' 


< 1 


III 


4 


3 




SA2For-' 
PO 7kXo 


I 36A2For-' 
PO ^ 7kXo 


-IOA2 


fo - abFoR'/''' 
where 

-205x2^''"' 
" ~ 7x(5)^/i» ' 

6 = (^)3/io 


< 1 


IV 


6 


2 


A2 


2A2For— 
PO kXo 


, 12A2For— 
PO ^ kXo 


28A2 
Xor^'^ 


fo - a5Foi?5/6 
where 

_ 3x2''/3 
" 75/6 ) 


< 1 


V 


6 


5 


A2 


2A2F„r-' 
PO 7kXo 


-Po+'-^ 


28A2 


fo - abFoR^^' 
where 

a - 

5^(A^)5/12 


< 1 


VI 


8 


5 


A2 


2A2For — 
P^ K.X0 


-P^^'^^l^ 


46A2 
Xor^k 


fo + abFo^ 
where 

123 

" 2-V3 X (23)2/3 ' 

b = {^Y" 


< 1 



The solutions corresponding to I2 are shown in Table 4. In this table, the increment of 
radius monotonically decreases pressure and density. Moreover, sound speed is subluminal. 
Thus we can conclude that the most acceptable solution of a cylindrical star in f{R) gravity 
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is generally described as: 



A{r) = A2r\ 

p{r) = Po+l>2r'+"-™-^ 

p{r) = -po + P2r'+"-"^-^ 

R{r) = i^2r'-™-^ R2 = ^i2f + iQ-m)l + 4m), 

f{R) = /o + a6Foi?^+^ k= " 



l-m-2 



If 3n 

a = - J (2/2 + (6-m)/ + 4m)-— (/-2m-4)-2n(n-l) 



(10/ + 6n-6m- 4){2n^ - 2n 
I + n — m — 2 



h = {^f^ , l = ^(m-2n- V32 + 16m + + 64n + 12mn - 28n^^ . 

4 Summary 

This paper is devoted to study the static cylindrically symmetric solutions in metric f{R) 
gravity for perfect fluid. In this theory, the field equations are highly non-linear and com- 
plicated which cannot be handled analytically without taking some assumptions. For this 
purpose, we confine ourselves to only those solutions which are closed and are found ana- 
lytically without use of any numerical data, perturbation or approximation method. The 
corresponding form of f{R) is reconstructed which does not depend algebraically on R{r) 
and /(r). Physical acceptability criteria is applied to check the physically acceptable solu- 
tions. We can summarize the results as follows: 

• First type of assumption provides some solutions given in Table 1 but none of them 
satisfies physical acceptability measures. Thus no new solution is obtained. 

• Second type is based on constant scalar curvature condition which gives some new 
solutions in Table 2. It is noted that some of them are physically acceptable but do 
not fulfill the imposed conditions. 

• For the third family of solutions, we formulate two tables corresponding to the pa- 
rameters li and l2- We obtain one solution in Table 3 which satisfies the physically 
acceptable criterion, whereas all the solutions in Table 4 fulfill the conditions for a 
physically acceptable solution. The later gives the static cylindrically symmetric in- 
terior solutions in the general form describing properties of a cylindrical star for f{R) 
gravity. 

Since the exact solutions of f{R) cosmological models are used to explore dark energy 
and dust matter phases. Thus, these solutions may provide a pointer towards the unknown 
nature of dark energy and dark matter. These solutions might be useful at some stage 
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to overcome the theoretical difficulties in the context of cosmology and astrophysics. The 
applicability of the solutions could be tested by comparing with cosmology and some other 
constraints. It would be interesting to investigate solutions for non-static spacetimes and 
also with different types of fluid. 
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